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基于 几何 形态 群 组 特征 的 R 波 检测 算法 


Xj Bs E M5, cH. Ax 维 ***， 蒋 建物 “ 
(郑州 大 学 a 软 件 与 应 用 科技 学 院 ;b. 产 业 技术 研究 院 ; c. 互 联网 医疗 与 健康 服务 河南 省 协同 创新 中 心 ， 郑 州 450000) 


摘 要 : 准确 检测 心 电 R 波 是 辅助 诊断 心脏 疾病 的 关键 。 针 对 OR 波 检测 问题 ， 提 出 了 基于 几何 形态 群 组 特征 的 及 
波 检 测算 法 。 算 法 首先 用 滤波 器 组 对 心 电 信 号 进行 预 处 理 ， 然 后 通过 整体 与 局 部 兼顾 的 信号 处 理 方式 计算 几何 形态 
特征 值 ， 并 使 用 群 组 中 三 个 几何 形态 特征 值 检 测 及 波 ， 最 后 采用 MITBIH 标准 心律 失常 数据 库 验 证 算法 。 实 验 结果 
表明 ， 算 法 对 及 波 检测 的 错 检 率 为 1.07%， 灵 敏 度 为 99.38%， 精 确 率 为 99.61%， 准 确 率 为 98.99%。 

关键 词 : 滤波 器 组 ; 几何 形态 ; 群 组 特征 ; 有 R 波 检测 
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R wave detection algorithm based on geometric morphology group characteristics 


Liu Wei* ^ °, Wang Lin® °, Li Runchuan" °, She Wei? ® °, Jiang Jianyu*' 
(a. Software & Applied Science & Technology Institute, b. Industrial Technology Research Institute, c. Internet medical & 
health service Henan Collaborative Innovation Center, Zhengzhou University, Zhengzhou 450000, China) 


Abstract: The accurate detection of ECG R wave is the key to assist in the diagnosis of heart disease. For the R wave 
detection problem, this paper proposes an R-wave detection algorithm based on geometric group features. Firstly, the 
algorithm preprocesses the ECG signal with the filter bank. Then, it calculates geometric eigenvalues by means of signal 
processing giving consideration to both overall and local, and uses the three geometric eigenvalues in the group to detect the 
R wave. Finally, it adopts MIT-BIH Standard Arrhythmia Database to test and verify this algorithm. The experimental 
results show that the error detection rate of the algorithm for R wave detection is 1.0796, the sensitivity is 99.3896, the 
accuracy rate is 99.6146, and the precision rate is 98.99%. 


Key words: filter bank; geometric morphology; group characteristics; R wave detection 


量 大 , 处 理 大 量 数据 时 实时 性 差 。 Chen 等 人 00 提 出 的 基于 自 

适应 阔 值 的 心 电 尺 波 检 测算 法 ,通过 设置 阔 值 并 且 不 断 检 查 
疾病 也 逐 每 个 心动 周期 的 方式 来 提高 算法 的 准确 率 ， 实 验 结果 显示 该 
体现 ,能 。 算法 提高 了 R 波 检测 的 准确 率 。 动 态 自 适应 闵 值 法 虽然 提高 


0 引言 


心血 管 疾病 严重 威胁 人 类 健康 ， 在 我 国 心血 管 
渐 成 为 高 发 病 。 心 电 图 (ECG) 是 心脏 活动 的 综合 


够 反映 人 体 心 脏 的 健康 状况 , 革 些 病变 可 以 通过 ECG 波形 的 ”了 检测 的 准确 率 ， 但 是 每 次 检测 都 要 对 阔 值 进行 更 新 ， 增 加 
异常 体现 出 来 02, 目前 心电图 广泛 应 用 在 医学 上 辅助 医生 鉴 。 了 算法 的 复杂 度 。 
别 与 分 析 各 种 心脏 疾病 。 心 电波 形 的 检测 在 心电图 分 析 上 起 在 以 上 方法 的 基础 上 ， 许 多 学 者 通过 改进 常用 算法 与 结 
通过 心 电 波形 的 检测 可 以 得 到 准确 的 心 。 合 不 同 的 算法 提出 了 一 些 新 的 方法 来 进行 R atem, 但 
重要 指标 。 ECG 检测 能 够 及 早 ”是 很 少 能 够 兼顾 到 算法 的 准确 率 、 普 适 性 和 简便 性 。 本 文通 
的 发 现 心脏 病 ,通过 预防 和 治疗 可 以 降低 心脏 病人 的 死亡 率 。 ”过 分 析 常 用 有 波 检测 方法 的 优 缺 点 ， 在 Yazdani 等 人 na 提出 
在 ECG 的 各 种 波形 中 R 波 是 心 电 信号 分 析 的 基础, 是 检测 心 。 的 自 适应 数学 形态 学 检测 算法 的 启发 下 提出 了 基于 几何 形态 
律 是 否 正常 以 及 诊断 心血 管 疾病 的 关键 。 群 组 特征 的 R 波 检测 算法 ,使 用 MITBIH 标准 心律 失常 数据 
差分 阅 值 法 是 最 早 应 用 在 R 波 检测 上 的 算法 , 也 是 R 波 “， 库 的 心 电 信号 数据 作 样 本 数据 进行 实验 ， 实 验 结果 表明 本 文 
检测 的 经 典 方法 ， 但 是 当 出 现 强 干扰 时 算法 的 检测 效果 会 有 提出 的 算法 具有 较 高 的 灵敏 度 和 准确 率 而 且 普 适 性 好 ， 运 算 
明显 的 下 降 。 随 着 计算 机 技术 的 发 展 与 科研 人 员 对 心 电 信号。 简单。 
检测 研究 力度 的 不 断 增 大 ， 各 种 新 型 的 信号 处 理 算法 也 被 应 
用 在 心电图 的 R 波 检测 中 。 常 用 算法 有 小 波 变换 法 Ba, 动 ” 1 ”数据 预 处 理 
d EDAM, 神经 网 络 U 等 方法 。 例 如 Razavit 等 人 人 体内 心 电 信号 频率 低 ， 容 易 受到 各 种 因素 的 干扰 。 因 
提出 基于 小 波 变换 的 方法 检测 R 波 ， 先 对 心 电 信号 进行 预 此， 通常 采集 的 心 电 信号 会 掺 杂 着 各 种 咒 声 。 相 同 的 信号 源 
处 理 滤 除 噪声 ， 然 后 对 信号 进行 连续 小 波 变换 ， 同 时 采用 阔 “产生 的 心 电 噪声 对 心 电 信号 的 影响 具有 相似 性 ， 通 常 可 以 把 
值 法 来 检测 QRS 波 群 的 R 波 ， 实 验 结果 表明 该 算法 有 较 高 。 这 些 噪声 归 为 一 类 ， 然 后 选取 合适 的 方法 进行 滤波 。 典 型 的 
的 准确 性 。 基于 小 波 变 换 的 检测 算法 虽然 检测 率 高 ,但是 运算 “噪声 有 肌 电 干扰 、 基 线 漂移 ， 肌 电 干扰 属于 高 频 干 扰 品 声 
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基线 漂移 属于 低频 噪声 。 通 常 进行 R 波 检测 时 都 需要 先进 行 。 1.2 中 值 滤波 器 

滤波 ， 滤 除 心 电 信号 中 的 噪声 对 于 R 波 检测 非常 关键 。 对 波 中 值 滤波 是 以 排序 统计 理论 为 基础 的 有 效 抑制 噪声 的 信 

波 后 的 心 电 信号 再 进行 R 波 检测 ， 能 够 减少 错 检 的 发 生 ， 提 ”号 处 理 技术 。 一 维 信号 处 理 最 先 应 用 中 值 滤波 技术 ， 随 后 二 

高 波形 检测 的 准确 率 。 维 信号 处 理 也 逐渐 开始 应 用 。 中 值 滤波 的 原理 是 把 数字 图 像 
常用 的 高 频 滤波 器 有 FIR 低 通 滤波 器 和 带 通 滤波 器 ,FIR ”或 数字 序列 中 某 一 点 处 的 值 用 该 点 的 一 个 邻 域 中 各 点 值 的 


中 
低 通 滤波 器 能 够 在 保证 任意 幅 频 特性 的 同时 具有 严格 的 线性 。”” 值 代替 。 设 一 维 序列 为 *,%, 为 … 为 ， 把 这 些 元 素 按 从 小 到 大 


相 频 特性 ， 又 因 其 具有 有 限 长 的 单位 抽样 响应 ， 因 此 滤波 器 。 排列 进行 中 值 滤波 ， 如 果 元 素 总 数 为 偶数 ， 则 中 值 为 中 间 
有 稳定 性 在 图 像 处 理 领域 中 应 用 广泛 。 带 通 滤波 器 的 原理 元素 的 均值 ， 元 素 总 数 为 奇数 ， 中 值 为 中 间 的 元 素 值 P0， 公 
是 设置 上 限 、 下 限 两 个 频率 界限 , 输入 信号 的 频率 处 于 上 限 、 — 式 如 下 : 

下 限 之 间 才 能 够 通过 滤波 器 ， 超 出 或 者 低 于 设置 的 值 都 将 会 i o 

被 过 滤 掉 。 带 通 滤波 器 整个 频率 通 带 内 都 能 最 大 限度 的 让 信 a aa o ee Pru] tim 

号 通过 ， 而 其 阻 带 对 信号 有 衰减 和 抑制 作用 。 基 线 漂移 噪声 "P 

的 滤 除 可 以 使 用 中 值 滤波 器 nb9 和 多 项 式 拟 合 滤波 器 。 中 值 滤 = 

波 器 运算 速度 快 并 且 设计 简单 ， 对 于 使 用 整数 表示 的 图 像 ， 记 为 y, e medix,] (8) 
因为 从 窗口 到 窗口 更 新 直方 图 简单 而 且 寻找 直方 图 的 中 值 不 y, — HERE BU somos 的 中 值 。 把 中 值 滤波 的 概念 
繁琐 ， 所 以 中 值 滤波 器 更 加 有 效 。 多 项 式 拟 合法 先 对 信号 进 ”推广 到 二 维 ， 选 用 特定 形式 的 二 维 窗口 对 二 维 序列 进行 中 什 


行 多 项 式 拟 合 ， 然 后 用 原始 信号 减 去 拟 合 信号 ， 得 到 去 除 基 ， 滤波， 二 维 信号 的 中 值 滤波 公式 如 下 ; 

线 的 信号 ， 但 是 多 项 式 拟 合 在 处 理 较 大 的 漂移 时 效果 不 佳 。 六 =medsfX 让 4 为 滤波 窗口 (9) 
通过 对 比 各 类 滤波 器 的 优 缺 点 ， 本 文 最 终 采 用 KEJE "—" 

通过 对 比 各 类 滤波 器 缺 本 文 最 终 采 用 FIR 低 通 滤波 2 QGMGC R 算 法 


器 先 对 心 电 信 号 的 肌 电 干扰 噪声 进行 滤 除 ， 然 后 用 中 值 滤波 

器 对 心 电 信号 再 次 滤波 ， 消 除 基线 漂移 噪声 ， 由 FIR 低 通 滤 本 文 提 出 的 基于 几何 形态 群 组 特征 的 R 波 检测 算法 

波 器 和 中 值 滤波 器 构成 滤波 器 组 进行 信号 预 处理 。 GMGC R 的 创新 点 为 对 信号 的 处 理 做 到 了 整体 与 局 部 的 兼 

1.4 FIR 滤波 器 厌 、 自 适应 性 好 。 整 体 与 局 部 结合 的 信号 处 理 方式 提升 了 检 
在 数字 信号 处 理 中 FIR 滤波 器 应 用 广泛 ， 主 要 功能 是 把 。 测 的 准确 率 ， 在 对 信号 求 几何 特征 值 了 = 时 ， 某 一 


干扰 信号 减弱 或 滤 除 得 到 所 需 的 信号 。 本 文 设计 的 基于 窗 函 。 mH Voss v 如 果 远 远 高 于 其 相 邻 点 的 了 = 
数 的 FIR 低 通 滤波 器 就 是 有 选择 性 的 减弱 或 滤 除 一 些 特定 频 。” 则 会 影响 算法 中 几何 形态 的 判定 进而 影响 算法 的 准确 率 ， 本 
率 范围 内 的 信和 号。 从 时 域 出 发 利用 窗 函 数 法 设计 FIR 滤波 器 。 ”文通 过 计算 连续 四 点 V=f,%,,…sv} HIE V merae] AE 
通过 选取 合适 的 窗 函 数 ， 把 无 限 长 的 理想 (0) 截 为 有 限 长 的 降低 茶 一 异常 点 对 算法 性 能 的 影响 。 在 几何 形态 的 趋势 特征 
AD ， 用 xD 逼近 为 (m) ， 使 频率 响应 HC) 和 所 求 的 理想 频 — 值 计算 中 ， 先 找 出 几何 形态 上 升 或 下 降 趋势 区 间 内 
率 响 应 Hu(e*) 接 近 。 其 中 窗 函 数 的 选取 对 于 设计 FIR 滤波 器 V S WiVa Va) RIEA AE V S Piva} ， 然 后 求 出 此 区 间 内 
是 很 重要 的 ， 需 要 通过 所 要 设计 的 滤波 器 的 性 能 指标 对 窗 函 V = [MYV HIJE V S Dovo}, died Voies) 
数 的 类 型 以 及 窗口 的 长 度 进行 确定 ,最 后 进行 整体 性 能 检测 。 I Vm Miva v) 的 均值 作为 几何 形态 的 趋势 特征 值 , 这 样 运 
数字 滤波 器 传输 函数 为 He"), h 是 与 理想 频 响 对 应 算 的 目的 也 是 为 了 消减 异常 点 对 算法 性 能 的 影响 ， 经 过 两 次 


的 理想 单位 抽样 响应 序列 ， 公 式 如 下 : 运算 ， 消减 了 异常 点 对 算法 性 能 的 影响 ， 提 升 了 R 波 检测 的 
"E "m 准确 率 。 
H,(e")- 2 h e” (1) E ee LLL ERES 
=o E3388 2J] 8s A EMRE In], ARAE E Y SLA 


TECTUM o POR, RESUH UI RIER BER, 
iiit E 的 是 为 了 减少 不 必要 的 几何 形态 判别 提高 算法 的 效率 ， 而 不 
[Ls A, (n) 通常 是 无 限 长 的 非 因果 序列 , 而 本 文 需要 的 为 。 需要 在 进行 R 波 检测 时 再 动态 的 改变 阔 值 。 然 后 由 相 邻 且 连 
有 限 长 的 低 通 滤波 器 ， 为 了 达到 设计 目的 ， 可 以 采用 对 ” 续 的 三 个 几何 形态 组 成 一 个 群 组 ， 通 过 群 组 中 几何 形态 的 特 


E No 


局 (进行 近似 逼近 的 方法 。 常 用 方法 是 用 窗 函 数 法 进行 MET R 波 检测 。 算 法 在 确定 R 波 位 置 的 思想 上 与 医生 从 

FIR 低 通 滤波 器 设计 Po,， 选取 适当 的 窗 函 数 截取 访 (， 公 式 。 视觉 角度 直观 的 从 心电图 中 对 R 波 进行 定位 一 致 。 现实 生活 

如 下 中 医生 通过 心电图 可 以 直观 的 找 出 幅 值 最 大 的 波形 ， 然 后 通 

h(n) 2 (n)wn) G) — 过 判断 此 波形 的 幅 值 是 否 高 于 其 前 后 两 个 相 邻 波形 的 幅 值 就 

ER GO 中 ，w(n) 为 偶 对 称 且 有 限 长 的 窗 ， 本 文 设计 。 可 以 确定 了 波 的 位 置 。 本 文 提出 的 算法 中 由 连续 三 个 几何 形 
时 选用 的 窗 函 数 为 汉 明 窗 ; 态 组 成 一 个 群 组 ， 通 过 群 组 中 几何 形态 的 特征 值 

vin) [ost-osscos( 2E JI m " wu ee e 

4 特征 值 来 蔡 代 心电图 中 直观 的 幅 值 ， 如 果 位 于 群 组 中 间 的 

其 幅度 响应 为 几何 形态 特征 值 同 时 大 于 其 前 后 两 个 几何 形态 的 特征 值 则 可 


EDNA] a] (5) 以 确定 中 间 的 几何 形态 为 候选 波 。 算 法 进行 R 波 检测 时 不 
Ne 需要 每 次 动态 的 改变 闵 值 来 提高 检测 的 准确 率 ， 降 低 了 算法 
w) 对 称 性 以 及 卷 积 定理 得 如 下 公式 : 的 复杂 度 ， 提 高 了 自 适 应 司 时 算法 运用 医生 直观 检测 R 
波 的 思想 以 群 组 中 几何 形态 的 特征 值 代 蔡 实 际 幅 值 来 检测 R 
波 而 不 是 以 心 电 信号 中 芯 点 对 R 波 进行 检测 ， 对 异常 波 
窗 函 数 法 是 应 用 窗 函 数 取 有 限 长 脉冲 响应 ， 以 4(m) 来 表 形 也 可 以 做 到 有 效 检测 ， 因 此 本 文 算法 具有 普 适 性 。 算 法 描 
ÀE h(n) 的 方法 ， 在 FIR 滤波 器 设计 中 应 用 广泛 。 述 如 下 所 示 : 
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4. while Dz(Z 
end 
5. while Vð, k-Vi-V , 
k, 2(,,7V)/3 
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S«—-SUK, , end 
6. while Cx , Wi 
同 理 得 Kyo, ， 
K«-KUK, , end 
y while CØ, 
if K «Ki, Ku, < Kao 
8. while Y z, 
if 
R«-RUY,, end 
9. return R 
本 文 算法 中 输入 的 是 心 电 


是 RR 波 位 置 的 集合 R。 


首先 把 集合 V H 


对 值 运算 ， 把 运算 后 上 
值 做 一 阶 差分 运算 , d 


U,, 的 正 负 ， 确定 几 n 


的 结 
巴结 果 存 入 数组 UV。 然后 依次 判断 
峰值 点 56, ,结束 点 6。 


ÜH NN 
PA 
ES 


的 起 始点 6,， 


直至 取 完 Uv 中 的 值 。 
形态 的 起 始点 6,， 
态 峰值 点 6,， 接 着 判断 UV, M Un, ÆU FEH 
则 该 点 为 这 个 几何 形态 的 结束 点 G. ， 同 时 也 是 后 继 几 何 
ESSE Ga 6,. G 存 入 数组 D n. Xt 
中 相 邻 且 连 续 的 起 始 、 峰 值 、 结 束 三 点 确定 一 个 完整 的 几何 


的 起 始点 ， 将 几何 


E 


fr 


数组 S H 
IN o 


Ph， 每 个 几何 


fiU, JV FE BU 大 于 零 则 该 点 为 几 
U, K-F38 HU. 小 于 零 则 该 点 为 


8*9 (ECG) 集合 VY ， 输 出 的 
FP 的 所 有 数值 进行 取 绝 
果 存 入 集合 v 中 。 对 集合 V 中 所 有 


IN 


U, fl 


IL 何 
U,, X TY 
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Ù wW 


计算 几何 


形态 


IZ RS s 


几何 特 4 


的 ] 


Efü k, , 
几何 特征 最 大 ff 
匀 值 Ka ， 最 后 计算 几何 
上 述 方式 计算 得 到 下 降 趋 势 特征 值 Cow o 
PARIRI Ke MH 
数组 C 与 数组 K 下 标 一 一 对 


ESSO 


并 从 数组 S FAH 


与 


取 每 个 点 后 连续 三 点 的 值 


形态 C ， 把 找 出 的 几何 形态 存 入 数组 C。 设 置 一 个 最 低 有 效 
装 值 工 即 电压 为 0.01 mV， 几 何 形态 特征 值 为 0.001， 滤 除数 
组 C 中 无 效 几 何 形态 ， 得 到 滤 除 无 效 几何 形态 后 的 数组 C， 
同时 把 集合 V 中 对 应 无 效 几何 形态 的 值 也 删除 得 到 删除 值 


后 的 集合 。 以 一 阶 微分 为 基础 ， 对 集合 Y 中 的 每 点 求 几何 
特征 值 Ks ， 即 从 集合 Y 中 依次 ] 
分 别 求 出 Ks 值 , 然后 取 这 四 个 Ki, 值 的 平均 值 
多 态 都 与 S 中 的 一 些 几 何 特征 值 相 对 


TER EA 


C BUB 


然后 从 这 些 中 找 出 起 始点 至 


EE Ka, MZH C 中 依次 取出 几何 
的 几何 形态 相对 应 的 所 有 


峰值 点 区 间 的 


直 Ku ， 并 计算 出 这 


P As Jy 0.6, JE Ke 值 存 入 数组 


M. 


区 间 中 


多 态 的 上 升 趋势 特征 值 Ko o 


所 有 几何 特 人 


K, FI ku 算出 


选 出 候选 R W, 依次 比 对 数组 C 中 的 三 个 连续 几何 形态 


的 Ke 值 ， 把 每 个 几何 


ESTO 


应 的 K şs 值 记 为 Kj. Ka. Kp, 


如 果 K; 小 于 Ks H Ki 大 于 Ks 时 ， 将 Ci 记 为 候选 R 波 ， 
可 形态 的 所 有 信息 存 入 数组 了 中 ， 信 息 中 包含 几何 形 
值 。 确 定 有 效 R 波 ， 依 次 比 对 数组 了 中 的 三 个 连续 


EH K. 值 ， 


将 此 几 
SM K, 


JL fij 


ade. 
TP FP 


F 1/2 WY HY, 大 于 1/2 的 %; 时 ， 则 可 以 判定 
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SEN = x 1009 K 
把 sev — T aos 值 记 为 
TP 
= 一 一 一 一 一 x100% . Yo, WE Ya K 
TP+FP+EN 6s Ya, WA Ka 


ES 
YE TAUR 


波 ， 将 其 存 入 数组 R， 输 出 即 为 R 波 位 置 的 集合 。 


续 处 理 得 出 最 终结 果 ， 而 且 算 法 中 的 每 个 步骤 主要 是 通过 


本 算法 采 上 


只 处 理 其 对 应 阶段 的 问题 


层 循序 遍历 来 处 理 前 一 结 


合 存储 数据 ， 
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3. 


电信 号 用 滤波 器 组 进行 预 
进行 R 波 检 测 。 RE V 
对 值 运算 
果 找 出 几 


天 


流水 线 式 的 处 理 模 式 ， 算 法 中 的 每 一 个 步 又 
得 出 中 间 结 果 ， 后 续 步 又 依次 连 


果 数 据 集 。 文 章 中 算法 采 | 
此 综合 上 述 分 析 本 文 提出 的 GMGC_R 算法 


CES 


数组 


有 较 低 的 运算 复杂 度 ， 


实验 结果 


1 实验 流程 


其 整 体 复杂 度 为 O(n)。 


本 文 按 图 1 所 示 流 程 进行 实验 验证 算法 。 原 始 采集 的 心 


中 的 


处 
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形态 的 上 升 趋势 特征 值 。 同 


公式 计算 几何 形态 
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心 电 数据 预 处 理 
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1 实验 流程 


Fig. 1 


Experimental flow chart 
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3.2 结果 与 讨论 

心 电 信 号 检测 算法 的 性 能 需要 权威 的 心 电 检 测 标 准 来 对 
其 进行 验证 ， 新 的 算法 只 有 经 过 权威 心 电 数 据 库 的 验证 后 才 
能 够 得 到 人 们 的 认可 。 通常 用 MIT-BIH 数据 库 对 提出 的 检测 
算法 进行 验证 。MITBIH 库 记 录 了 48 个 病例 ,每 个 病例 数据 
长 30 分 钟 ， 总 计 约 有 116000 多 个 心 拍 ， 其 中 既 有 正常 心 拍 
出 有 有 异常 心 拍 ， 可 以 对 及 波 检测 算法 进行 有 效 的 验证 C2。 
本 文 所 有 实验 都 在 相同 实验 环境 下 进行 ， 其 中 CPU 为 
Intel®Core™i7-7700@3.60 GHz, 内 存 为 16 GB， 操作 系 统 为 
64 位 Windows10。 在 MATLABR20I5b 中 用 滤波 器 组 对 心 电 
数据 进行 预 处 理 , 然后 用 GMGC. R 算法 进行 R 波 检测 。 MIT 
心 拍 定位 标注 包括 N(1)、L(2)、R(3)、a(4)、V(5)、F(6)、J(7)、 
A(8)、 S(9), ECOS j(1D. /(12)、 Q(13), |(16)、!(31)、 e(34). 
n(35)、f(38)， 检 测 结果 如 表 1 所 示 。 表 1 中 的 检测 指标 有 四 
项 ， 错 检 率 记 为 DER， 灵敏 度 记 为 SEN， 精确 率 记 为 +P， 准 
确 率 记 为 Acc， 公 式 如 下 所 示 : 


T 


DER = ETEN X100% (10) 
TP 

SEN = — E — x100% (11) 
TP FN 

+P =P x100% 2) 
TP+ FP 

Acc=— P — 0096 (13) 

TP FP FN 


n 


€ 1 MITBIH 心率 失常 数据 库 R 波 的 检测 部 分 结果 统计 
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确 率 、 准 确 率 为 100%， 主 要 原因 是 这 类 ECG 信和 号 自身 携带 
的 噪声 较 少 ， 前 期 滤波 处 理 彻 底 滤 除了 噪声 ， 如 100 号 、103 
号 和 109 号 等 。 通 过 表 1 中 实验 结果 可 以 看 出 ， 除 了 207 号 
文件 因 噪声 过 多 造成 实验 结果 不 理想 外 ， 该 算法 每 一 个 文件 
的 灵敏 度 为 96.36% 以上， 精确 率 为 95.30% 以 上 ， 准 确 率 为 
94.43% 以 上 。 所 有 文件 的 错 检 率 为 1.07%, 灵敏 度 为 99.38%, 
精确 率 为 99.61%， 准 确 率 为 98.99%。 通 过 每 一 个 文件 各 项 
指标 的 结果 和 总 体 文件 的 检测 结果 可 以 看 出 本 文 算法 不 但 能 
够 准确 检测 出 正常 波形 的 R 波 ， 对 异常 波形 的 R 

波 检测 也 具有 较 好 的 效果 ， 说 明 该 算法 具有 较 高 的 普 适 
性 。 将 本 文 R 波 检测 算法 的 实验 结果 和 文献 中 的 检测 算法 结 
果 进 行 对 比 ， 不 同 算法 实验 结果 对 比如 表 2 所 示 。 

表 2 算法 结果 对 比 表 


Table 2 Algorithm result comparison table 


算法 SEN +P Acc 
文献 [8] 99.3990 99.4996 98.89% 
文献 [12] 98.4696 98.6296 97.12% 
文献 [13] 98.6896 97.24% 96.03% 
文献 [16] 96.2896 99.71% 97.65% 
文献 [17] 98.89% 99.40% 98.37% 
This paper 99.38% 99.61% 98.99% 


Table Statistical results of R wave detection in MIT-BIH arrhythmia 
database 
FileName DER SEN +P Acc 
100 0 100.0096 100.0096 100.0096 
101 0.5996 99.6396 99.79% 99.41% 
102 0.74% 99.36% 99.91% 99.27% 
103 0 100.00% 100.00% 100.00% 
106 0.25% 99.85% 99.90% 99.75% 
107 0.23% 99.86% 99.91% 99.77% 
109 0 100.0096 100.0096 100.0096 
113 0 100.0096 100.0096 100.0096 
114 0.0596 99.95% 100.00% 99.95% 
117 0 100% 100.00% 100.00% 
119 0.05% 100.00% 99.95% 99.95% 
124 0.12% 100.00% 99.88% 99.88% 
202 0.14% 99.86% 100.00% 99.86% 
205 0.15% 99.85% 100.00% 99.85% 
207 14.92% 87.32% 99.61% 87.02% 
212 0.25% 100.00% 99.75% 99.75% 
213 0.18% 100.00% 99.82% 99.82% 
217 0.41% 99.82% 99.77% 99.59% 
220 0 100.0096 100.0096 100.0096 
231 0 100.0096 100.0096 100.0096 
232 0.0696 100.0096 99.9496 99.94% 
Total avg 1.07% 99.38% 99.61% 98.99% 
本 文 用 四 个 指标 对 算法 性 能 进行 衡量 ， 由 表 1 中 可 以 看 


出 算法 的 错 检 率 为 1.07% ， 有 灵敏度 为 99.38%， 精 确 率 为 
99.6190, 准确 率 为 98.99%。 在 实验 结果 中 出 现 了 207 号 文件 
灵敏 度 为 87.32%， 准 确 率 为 87.02% 的 现象 ， 主 要 原因 是 这 
类 心 电 信号 自身 带 有 较 强 的 噪声 ， 经 过 滤波 器 组 滤波 后 ， 没 
有 彻底 滤 除 噪声 。 其 中 部 分 心 电 信 号 R 波 检测 的 灵敏 度 、 精 


从 表 2 中 可 以 看 出 与 其 他 文献 中 提出 的 及 波 检测 算法 相 
比 ， 本 文 算法 具有 较 好 的 灵敏 度 、 精 确 率 和 准确 率 。 本 文 算 
法 核心 部 分 是 以 整体 与 局 部 兼顾 的 信号 处 理 方式 计算 几何 形 
态 特征 值 ， 算 法 处 理 信 号 时 兼顾 整体 与 局 部 的 思想 ， 既 可 以 
规避 因 某 一 点 几何 特征 值 偏 大 对 实验 结果 造成 影响 又 具有 滤 
波 功能 。 兼 顾 整 体 与 局 部 的 运算 消减 了 异常 点 对 算法 性 能 的 
影响 ， 提 升 了 R 波 检测 的 准确 率 。 算 法 的 另 一 创新 点 为 用 几 
何 形态 群 组 特征 值 确定 R 波 , 通过 比 对 群 组 中 三 个 几何 形态 
特征 值 检 测 出 R 波 , 极 大 的 提高 了 本 算法 的 灵敏 度 和 精确 率 。 
综合 以 上 分 析 ， 本 文 提 出 的 GMGC_R 算法 在 灵敏 度 、 精 确 
率 、 准 确 率 、 普 适 性 方面 有 明显 的 优势 而 且 运 算 复杂 度 低 ， 
适用 于 各 种 心 电 信号 中 的 R 波 检测 。 


4 ARA 


本 文 提出 了 基于 几何 形态 群 组 特征 的 R 波 检测 算法 ， 首 
先 采 用 滤波 器 组 对 心 电 信 号 进行 预 处 理 为 准确 定位 R 波 英 定 
了 基础 。 然 后 在 确定 几何 形态 特征 值 的 运算 中 兼顾 了 整体 与 
局 部 信号 的 处 理 ， 不 但 降低 了 异常 信号 点 对 几何 形态 特征 值 
计算 的 影响 而 且 可 以 对 心 电 信 号 再 次 去 噪 ， 最 后 用 群 组 中 几 
可 形态 特征 值 对 波形 进行 检测 。 用 MIT-BIH 心律 失常 数据 库 
对 算法 进行 了 验证 ， 并 与 其 他 的 算法 实验 结果 做 了 对 比 ， 结 
果 表 明 本 文 提出 的 算法 有 较 高 的 灵敏 度 、 精 确 率 、 准 确 率 ， 


而 且 具 有 普 适 性 。 
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